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SUMMARY 

Exponential free-radical growth constants have been measured for 
formaldehyde-carbon monoxide-oxygen systems by monitoring the growth of 
oxygen atom concentration as manifested by CO flame-band emission. Data 
were obtained over the temperature range of 1200 to 2000 K. 

The data have been analyzed using a formaldehyde oxidation mechanism 
involving 12 elementary reaction steps. The computed growth constants are 
roughly in accord with experimental values, but are much more temperature 
dependent. The data have also been analyzed assuming formaldehyde is rapid- 
ly decomposed to carbon monoxide and hydrogen. Growth constants computed 
for the resulting carbon monoxide-hydrogen-oxygen mixtures have a tempera- 
ture dependence similar to experiment; however, for most mixtures, the com- 
puted growth constants were larger than experimental values. 

INTRODUCTION 

Formaldehyde is an important intermediate in the oxidation of hydrocar- 
bons to carbon dioxide and water. Thus, in order to obtain a complete set 
of chemical reactions and rates for modeling hydrocarbon combustion it is 
desirable to know the pathways by which formaldehyde reacts to form carbon 
monoxide, hydrogen, water, and carbon dioxide. 

The shock tube has proved useful for obtaining elementary reaction rates 
from growth constants in the H 2 -CO-O 2 (ref. 1) and CH 4 -CO-O 2 (ref. 2) 
systems. Consequently an analogous study of the CH 2 O-CO-O 2 system was 
initiated in hopes of extracting elementary reaction rates for formaldehyde 
oxidation. 


Formaldehyde Oxidation Scheme 

When a mixture containing carbon monoxide, oxygen, and a small amount of 
formaldehyde is subject to a temperature and pressure pulse in a shock tube, 
small concentrations of atoms and free radicals may first be formed by 
processes such as 





CO + 02“*^ CO 2 + 0 
CH 2 O + O 2 -*- CHO + HO 2 
CH 2 O + M— i^CHO + H + M 

These atom and free-radical concentrations may then grow exponentially via 
the branched-chain scheme 


^2 

H + O 2 — ^OH + 0 

(II) 

OH + C0_^C02 + H 

(V) 

kg 

0 + CH 2 O— §-CH0 + OH 

(VIII) 

^10 

CHO + M ^ CO + H 

(X) 

^11 

HO 2 + CO _i^C02 + OH 

(XI) 

k-4 

HO 2 + M H + O 2 + M 

Chainbranching is inhibited by the reactions 

(-IV) 

H + O 2 + M_4-H02 + M 

(IV) 

kfi 

OH + CH 2 O—S-CHO + H 2 O 

(VI) 

H + CH 2 O- 4 -CHO + H 2 

(VII) 

kg 

CHO + 02_^C0 + HO 2 

(IX) 

k-io 

H + CO + M -i^CHO + M 

(-X) 

ki 2 

0 + CO + M CO 2 + M 

(XII) 


The chain carriers H, OH, and 0 are more reactive than CHO, which in turn is 
more reactive than H02- Numerical integration of this kinetic scheme for 
the experimental conditions indicates that the rate of reaction (-IV) is 
greater than reaction (IV) and that there may be some reversal of reaction 
(V). Reactions have been numbered to be compatible with notation used in 
the study of the H 2 -O 2 -CO system (ref. 1). 

Theory shows (refs. 2 to 5) that in such chain branched systems when 
depletion of reactants is negligible and temperature and pressure are con- 
stant, the atom and radical concentrations increase propotional to exp 
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(xt) increases (except very early in the reaction). The growth constant 
X depends on the rate constants of the elementary chemical reactions and 
the concentrations of stable reactants - in this case CO, O 2 , and CH 2 O. 

This report presents experimental growth constants measured behind inci- 
dent shocks for a range of pressures, temperatures, and gas compositions. 
These experimental growth constants are compared with values computed from 
theory. Agreement between theory and experiment is only semi quantitative. 

Inasmuch as the analysis of the data is not definitive, the experimental 
conditions for each datum - gas time interval, temperature range, and pres- 
sure, are recorded for the benefit of others who may wish to reevaluate the 
data. 


EXPERIMENTAL ASPECTS 

Growth constants were obtained by measuring the blue CO flame band emis- 
sion behind incident shocks. The intensity of this radiation is proportion- 
al to the product of carbon monoxide and oxygen-atom concentrations 
(ref. 6), and since little or no CO is consumed, the light measures the in- 
crease of oxygen atom concentration with time. Details of the shock tube 
and associated optical and electronic equipment have been described else- 
where (ref. 7). 

Gas mixtures contained small amounts of CH^O with varying amounts of 
CO, © 2 , and COo diluted with argon. (Carbon dioxide was added to ensure 
vibrational relaxation of carbon monoxide.) Oxygen and argon were high- 
purity tank gases and were used without further purification. Carbon monox- 
ide was condensed at liquid nitrogen temperature; about a fourth of the con- 
densate was pumped off and discarded. Dry ice served as a convenient source 
of carbon dioxide. It was purified by subliming three quarters of a sample 
into a liquid nitrogen cooled trap. The first quarter of this trapped frac- 
tion was discarded and the middle half used for mixture preparation. 

Formaldehyde was prepared from trioxane, a cyclic trimer of formaldehyde 
which contains no combined water. This is important, since traces of water 
catalyze the polymerization of formaldehyde to paraformaldehyde. On the 
other hand, pure dry formaldehyde gas shows no visible signs of polymeriza- 
tion when stored in glass vessels at 80 to 100° C. 

Trioxane was depolymer ized by passing a nitrogen-trioxane mixture over 
an acid catalyst at 200 to 240° C. Formaldehyde was collected in a trap 
cooled by liquid nitrogen and purified by two bulb-to-bulb transfers in a 
glass-vacuum system. A 5 percent formaldehyde mixture in argon was prepared 
from the purified liquid and was stored in a glass bulb maintained at 
110° C. 

Test mixtures were prepared in a large glass vessel by introducing first 
the formaldehyde argon mixture and then adding an appropriate CO-O 2 -CO 2 
argon mixture from a high pressure storage tank. 

It has been shown (refs. 7 and 8) that boundary-layer effects must be 
considered in analyzing data obtained behind incident shocks. Conditions 
behind the shocks, in the region of the experimental measurements, were ob- 
tained from a computer program (ref. 9) which integrates the equations of 
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chemical change for a shocked-gas accounting for the effects of boundary 
layer buildup. The procedure has been described for both turbulent boundary 
layers (ref. 7) and laminar boundary layers (ref. 8). Data were obtained 
with both laminar and turbulent boundary layers in this report. 

Exponential growth constants were obtained from plots of logarithm of 
observed light intensity versus gas time. The relation between gas and 
laboratory times was obtained from the computer calculations. 

THEORETICAL CONSIDERATIONS 

The analytic solutions of the differential equations describing the ig- 
nition kinetics in chain-branched systems involving hydrogen, oxygen, and 
carbon monoxide have been discussed in detail previously (refs. 3 to 5) and 
the solution has been presented for the methane-carbon monoxide-oxygen sys- 
tem (ref. 2). The solution for the system of equations discussed in the 
introduction is quite analogous and will merely be outlined here. 

To obtain the analytic solution, the atom and radical concentrations H, 
OH, 0, HO 2 , and CHO are assumed to be small in comparison to the concen- 
trations of reactants CH 2 O, CO, and O 2 so that the reactant concentra- 
tions may be considered constant and reactions between chain carriers are 
unimportant. In addition, the temperature and pressure must be nearly con- 
stant over the time range of the experimental observations. 

Subject to these assumptions, the chemical kinetics are described by a 
system of simultaneous first-order linear differential equations, one for 
each chain carrier. The solution of these equations shows that the chain- 
carrier concentrations grow as exp (xt) (except very early in the 
reaction). The growth constant x is the positive root of a polynomial of 
a degree equal to the number of chain carriers. Thus, with H, OH, 0, HO 2 , 
and CHO as chain carriers, the following quintic equation is obtained; 


■(Vg 




-4 


10 




-(Vj ♦ 9) 


11 


0 

V. 




0 

V, 


-(v_ 4 ^Vii.X) 


-(Vg+v^Q-^^^ 


where = k2[02], e k^[02]CM], v_^ e 

V5 = l<5CC0], v_g = k_5[C02], vg = kg[CH20], vj = k7[CH20] 

^ '^9 ~ ^ g [^ 2 ^’ '*10 ^ '*10 " 


'*11 - 2nd v^2 = 
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The rows in this equation derive from the differential equations for the 
rates of formation of H, OH, 0, HO 2 , and CHO, respectively. The columns 
derive from the rates of formation or distruction of H, OH, 0, HO 2 and CHO 
in the elementary chemical reactions. 

This equation was used to calculate growth constants corresponding to 
the experimental conditions using rate constants taken from the literature. 
The rate constants and sources are listed in the appendix. Also, Equa- 
tion (1) was numerically differentiated to obtain the sensitivities of the 
calculated growth constants to the various reaction rates. 

RESULTS AND DISCUSSION 

The compositions of the four gas mixtures are given in table I. For 
mixture 1 there were slightly different compositions, designated as a, b, 
and c. Also shown are the sensitivities, a In x/a In for the various 

reaction rates. The range of sensitivities shown is for the temperature 

range of the data for each composition. On the basis of low sensitivities 

several reactions might be eliminated from the scheme reactions (IV), (-V), 

(-X), (XII), and perhaps (-IV). 

The experimental results are set forth in table II and Figures 1 to 4. 
The data for mixtures 1 and 2 were obtained at pressures (1 to 1.5 atm) 
where the boundary layers were turbulent. The data for mixtures 3 and 4 
were at lower pressures (0.12 to .25 atm) where boundary layers were 
laminar. The data for mixtures 1 and 2 range from about HOOK to ISOOK; the 
data for mixtures 3 and 4 cover a wider temperature range, from about 1200K 
to over 2000K. 

Because the most important reactions in the formaldehyde oxidation 
scheme are all bimolecular [only reactions (IV), (-X), and (XII) are not], 
the growth constant at a given temperature should be proportional to the 
pressure. Consequently, in figures 1 to 4 the experimental data are plotted 
as the logarithm of x/P versus reciprocal temperature. Also shown as 
solid lines are values computed from equation (1) using rate constants taken 
from the literature, and set forth in the appendix. 

The calculated growth constants are generally smaller than experimental 
values, except for mixtures 3 and 4 at the higher temperatures. Indeed, at 
the lowest temperatures, computed growth constants for mixtures 3 and 4 are 
about one seventh of the experimental values. 

Attempts were made to reconcile theory and experiment by changing the 
rates of several of the important reactions. From the sensitivities in 
tableil, the important reactions are (II), (V), (VI), (VII), (VIII), (IX), 
(X), and (XI). The rates of reactions (II) and (V) are well-established and 
are not candidates for tinkering. Values of kg were computed from mixture 
3 using equation (1) together with the experimental growth constants. Al- 
though a good correlation was obtained, when the rate const^ni was fitted to 
the Arrhenius equation the pre-exponential factor was 8X10^'^ which is 
three orders of magnitude too large for such a bimolecular abstraction 
reaction. In another attempt to fit the data k^o was computed from 
mixture 3; these calculations yielded rate constants which were larger at 
low temperatures than at high temperatures, which is absurd. In still 
another attempt to fix the data ky was computed from mixtures 1 and 2. 
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When these rate constants were fitted to the Arrhenius equation, the 
pre-exponential factor was about four orders of magnitude too large. 

The pattern which seems to emerge is the following: to fit the data, 

reactions that promote branching such as (II), (V), (VIII), (X), or (XI), 
must be assigned small or negative temperature dependences; or inhibiting 
reactions such as (VI), (VII), and (IX) must be assigned very strongly posi- 
tive temperature dependence. Neither alternative is reasonable. 

Perhaps the formaldehyde is rapidly decomposed to hydrogen and carbon 
monoxide so that the growth constants are really those of the H 2 -CO-O 2 
system. To explore this possibility, growth constants were calculated for 
the H 2 -CO-O 2 system, assuming the formaldehyde present decomposed in- 
stantaneously to hydrogen and carbon monoxide. Results are shown as dashed 
curves on figures 1 to, 4. For mixtures 3 and 4 (figs. 3 and 4), which span 
a wide range of temperature, growth constants computed for the H 2 -CO-O 2 
system show a temperature dependence similar to that of the experimental 
data; indeed, for mixture 3, calculated and experimental growth constants 
are in agreement. For mixture 4, calculated growth constants are perhaps 60 
percent to 90 percent larger than experiment. For mixture 1, experimental 
growth constants lie between values calculated for the CH 2 O-CO-O 2 and 
H 2 -CO-O 2 systems, while for mixture 2, calculated growth constants for 
both schemes are lower than experimental values. 

The analyses reported here suggest there may be substantial but incom- 
plete decomposition of formaldehyde to carbon monoxide and hydrogen prior to 
the exponential growth of chain carriers. To model this situation, numeri- 
cal chemical kinetic calculations would be required. It is not clear that 
such calculations would accomplish more than a rationalization of the exper- 
imental data; it is uncertain as to whether rates of elementary reactions 
could be established with confidence. 

CONCLUDING REMARKS 

Exponential free-radical growth constants have been measured for 
formaldehyde-carbon monoxide-oxygen systems by monitoring che growth of oxy- 
gen atom concentration as manifested by CO flame band emission. Data were 
obtained over the temperature range of 1200 to 2000 K. 

The data have been analyzed using a formaldehyde oxidation mechanism 
involving 12 elementary reaction steps. The computed growth constants are 
roughly in accord with experimental values, but are much more temperature 
dependent. The data have also been analyzed assuming formaldehyde is 
rapidly decomposed to carbon monoxide and hydrogen. Growth constants com- 
puted for the resulting carbon monoxide-hydrogen-oxygen mixtures have a tem- 
perature dependence similar to experiment; however, for most mixtures, the 
computed growth constants were larger than experimental values. 
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APPENDIX 


SPECIFIC REACTION RATE CONSTANTS 

This appendix lists individual reactions used in the calculations ana 
indicates values and sources of the assumed rate constants. Bimolecular 
rate constants are expressed in cm^mole“^sec”^ and the molecular rate 
constants are in cm°mole“^sec“^ . Temperatures are in degrees Kelvin 
and activation temperatures, E/R, are in degrees Kelvin. 

(II) H + O 2 OH + 0, k 2 = 1.25X1014 exp (-8203/T) 

This rate constant, from reference 1, was measured in the shock tube used in 
this investigation. 

(IV) H + O 2 + M + M, k 4 = l.SXlOlS exp (500/T) 

This rate constant (ref. 10) is for argon as the third body. Other third- 
body factors used (Ar = 1) were O 2 and CO = 1.3, CO 2 = 5, and CH 2 O = 1. 

(-IV) HO 2 + M-s^H + O 2 + M, k_4 = l.SXlolSexp (-25 260/T) 

This rate constant is obtained from the rate of reaction (IV) and the equi- 
librium constant for reaction (IV), from reference 11. Chaperon efficien- 
cies are the same as in reaction (IV). 

(V) CO + 0H-^C02 + H, k 5 = 10l2 exp (-1860/T) 

This rate constant (ref. 1) was obtained from the shock tube used in this 
investigation. 

(-V) CO 2 + H-i^CO + OH, k_5 = 1.24X1014 gxp (-13000/T) 

This rate constant is obtained from the rate of reaction (V) and the equi- 
librium constant for reaction (V) from reference 11. 

(VI) OH + CH20-5^H20 + CHO, ks = 7.8X10l2 
Recommendation of reference 12. 

(VII) H + CH20-i^H2 + CHO, ky = 1.35X10l3 exp (-1892/T) 

From reference 13. 

(VIII) 0 + CH20-^0H + CHO, kg = 1.77X1013 gxp (-1542/T) 

From reference 14. 

(IX) CHO + 02-5-CO + HO 2 , kg = 5X10l3 exp (-835/T) 

This rate constant is obtained by fitting the Arrhenius equation to the room 
temperature recommendation of reference 12 with the 1600K value of reference 
15. 

(X) CHO + M-i^CO + H + M, kio = 6.7X10l3 exp (-8741/T) 
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This rate constant is obtained from the rate of reaction (-X), and the equi- 
librium constant for reaction (X). The equilibrium constant was calculated 
from the thermodynamic data of reference 11, with an adjustment to account 
for the heat of formation of CHO recommended in reference 12. The rate 
shown is for M = Argon. Other chaperon efficiencies were taken to be those 
of reaction (IV). 

(-X) H + CO + M-j* 4:H0 + M, k_io = 10l4 exp (-755/T) 

This rate is 0.2 of that suggested in reference 16 with hydrogen as a third 
body and was assumed appropriate for M = Argon. Other chaperon efficiencies 
were assumed the same as reaction (IV). 

(XI) HO 2 + C0-i^C02 + OH, kii = 8 XIOI 3 exp (-11500/T) 

Recommendation of reference 16. 

(XII) 0 + CO + M-i»-C02 + M, ki 2 = 4X10l5 exp (-2300/T) 

Recommendation of reference 16. 
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TABLE I. - MIXTURE COMPOSITION AND GROWTH CONSTANT SENSITIVITIES 


Compositions, 
mole % 

Gas mixture 


1 

2 

3 

4 


a 

b 

c 




CHpO 

0.481 

0.530 

0.474 

0.117 

0.175 

0.170 

CO 

4.91 

4.96 

4.92 

5.36 

4.73 

3.32 


.520 

.525 

.474 

.567 

1.88 

15.30 

CO 2 

4.96 

5.01 

4.97 

5.42 

.964 

1.16 


Sensitivities, 
reaction i 

Gas mixture 

1 

2 

3 

4 

II 

1.12 to 1.05 

0.92 to 0.70 

0.96 to 0.35 

0.49 to 0.07 

IV 

-.02 to -0.01 

-.04 to -0.01 

-.01 to 0 

-.04 to 0 

-IV 

.05 to 0.06 

.01 to 0.02 

.01 to 0.05 

.02 to 0.19 

V 

.14 to 0.11 

.21 to 0.17 

.65 to 0.23 

.82 to 0.40 

-V 

-.03 to -0.05 

-.04 to -0.05 

0 to -0.02 

0 

VI 

-.10 to -0.02 

-.12 to -0.01 

-.63 to 0.07 

-.84 to -0.01 

VII 

-.89 to -0.62 

-.42 to -0.15 

-.89 to -0.04 

-.35 to -0.02 

VIII 

.04 to 0.07 

.19 to 0.24 

.03 to 0.32 

.05 to 0.33 

IX 

-.87 to -0.68 

-.43 to -0.22 

-.28 to -0.13 

-.08 to -0.12 

X 

.88 to 0.70 

.43 to 0.23 

.28 to 0.14 

.08 to 0.12 

-X 

0 

-.01 to 0 

0 

0 

XI 

.70 to -0.39 

.33 to 0.12 

.89 to 0.03 

.85 to 0.04 

XII 

-.02 to -.01 

-.06 to -0.03 

-.01 to 0 

-.01 to 0 
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TABLE II. - EXPERIMENTAL RESULTS 


Temper- 

ature, 

K 


Pressure, 

atm 


Growth 
constant, 
X ,sec“^ 


Temperature 

range, 

K 


Pressure 

range, 

atm 


Time 

range, 

usec 



1596 

1934 

1524 

1506 

1493 

1481 

1469 

1463 


1598 

1573 

1571 

1526 

1495 

1442 

1410 


1.1589 

1.2472 

1.1084 

1.1988 

1.3505 

1.3864 

1.3871 

1.4840 


1.0958 

1.1485 

1.1337 

1.1658 

1.2598 

1.3257 

1.4137 


26 . 62 x 103 

20.94 

17.12 

16.45 

18.28 

18.42 

16.93 

17.99 


31 . 99x10 ^ 

29.53 

30.50 

26.32 

23.99 

21.93 

23.38 


■ 


,576 - : 

514 - : 
1490 - 1557 
1480 - 1531 
1483 - 1503 
1460 - 1501 
1463 - 1474 
1457 - 1469 


Mixture 2 



■ 


.592 - 1603 
567 - 1578 
1563 - 1578 
1519 - 1533 
1490 - 1498 
1435 - 1448 
1405 - 1414 


Mixture 3 



150 - 240 
180 - 270 
150 - 260 
190 - 310 
230 - 305 
260 - 420 
380 - 510 


mmm 

0.2193 

11 . 000x103 

2038 - 2068 

0.204 - 0.216 

100 - 300 


.1472 

8.550 

2013 - 2041 

.145 - 0.149 

100 - 450 _ 

2022 

.2331 

11.250 

2010 - 2034 

.230 - 0.236 

75 - 300 

1832 

.2187 

12.250 

1823 - 1841 

.217 - 0.220 

150 - 400 

1774 

.1571 

7.425 

1757 - 1780 

.155 - 0.159 

200 - 700 

1728 

.1781 

7.533 

1719 - 1738 

.177 - 0.179 

350 - 675 

1669 

.2078 

8.200 

1658 - 1680 

.206 - 0.209 

300 - 675 

1618 

.1714 

8.220 

1609 - 1624 

.170 - 0.172 

525 - 850 

1611 

.1989 

7.967 

1602 - 1620 

.197 - 0.200 

275 - 700 

1611 

.1706 

5.580 

1600 - 1622 

.169 - 0.172 

550 - 1000 

1557 

.1401 

4.860 

1547 - 1567 

.138 - 0.142 

550 - 1150 

1351 

.2524 

3.125 

1343 - 1358 

.250 - 0.254 

1300 - 2300 

1342 

.1414 

3.224 

1338 - 1346 

.140 - 0.142 

1700 - 2750 

1297 

.1285 

2.305 

1293 - 1301 

.127 - 0.129 

2300 - 3900 

1296 

.1342 

2.515 

1291 - 1301 

.133 - 0.135 

2100 - 3900 

1266 

.1470 

2.320 

1263 - 1268 

.146 - 0.148 

2900 - 3950 

1258 

.1445 

2.940 

1253 - 1262 

.143 - 0.146 

1600 - 2750 

1211 

.1458 

1.870 

1208 - 1214 

.145 - 0.147 

3400 - 5200 
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Table II. - Concluded 


Temper- 

Pressure, 

Growth 

Temperature 

Pressure 

Time 

ature. 

atm 

constant, 

range. 

range. 

range. 

K 


x,sec"^ 

K 

atm 

usec 

Mixture 4 

2007 

0.1804 

9.202x10^ 

1989 - 2024 

0.178 - 0.183 

80 - 400 

1812 

.2260 

9.146 

1794 - 1830 

.223 - 0.229 

100 - 500 

1695 

.2750 

10.909 

1684 - 1706 

.273 - 0.277 

120 - 420 

1679 

.2147 

10.563 

1668 - 1690 

.212 - 0.217 

120 - 460 

1659 

.2255 

8.824 

1648 - 1670 

.223 - 0.228 

160 - 520 

1604 

.2444 

9.464 

1595 - 1612 

.242 - 0.247 

120 - 480 

1572 

.1917 

6.985 

1563 - 1580 

.189 - 0.194 

200 - 650 

1512 

.2148 

7.100 

1506 - 1517 

.213 - 0.216 

325 - 725 

1444 

.1992 

5.950 

1438 - 1449 

.197 - 0.201 

400 - 950 

1443 

.1782 

3.960 

1435 - 1450 
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